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1 Introduction 



The Large Hadron Collider (LHC) which will operate at an enormous center of mass 
energy (v^ = lATeV) offers to shed light on the mechanism of spontaneous symme- 
try breaking. It also promises to be a discovery machine and it is hoped that some 
signals of new physics beyond the standard model (SM) will be observed. Many 
exciting possibilities have been envisaged; most popular ones are supersymmetry, a 
symmetry which relates fermions to bosons, and the possibility of extra dimensions. 
In this paper we will consider the large extra dimension model by Arkani-Hamed, 
Dimopoulos and Dvali (ADD) |Ill|2l. There are many important discovery channels 
at the LHC such as 77, ZZ, W^W~ , jet production etc. These processes have already 
been studied in the context of anomalous triple gauge boson vertices ||3|1. In the SM the 
production of two Z bosons is suppressed as it begins at the order in the the electro- 
magnetic coupling and also because of the large ZZ production threshold. The two 
Z bosons can couple to Kaluza Klein (KK) gravitons, thus ZZ pairs can be produced 
through virtual graviton exchange at the leading order. These observations make ZZ 
production one of the important discovery channels. At LHC, Quantum Chromody- 
namics (QCD) plays an important role through higher order radiative corrections in 
reducing various theoretical uncertainties in the prediction of various important ob- 
servables to understand both SM as well as new physics scenarios. In this paper we 
will consider production of Z boson pairs at the LHC at next-to-leading order (NLO) 
accuracy in the strong coupling constant. 

Leading order studies for ZZ production in the SM can be found in [4j. Z pair 
with a large transverse momentum jet at LO was studied in JS]]. LO studies for ZZ 
production in the context of extra-dimension models have already been carried out 
in ©[Zll and limits on the parameters of the model have been placed. But, as we know 
a leading order (LO) computation is very sensitive to the choice of factorization scale 
(fip) which enters through the parton distribution functions at this order. To have a 
precise prediction of cross section and to have tighter constraints on the parameters 
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of the model it is necessary to carry out a full NLO calculation in the strong coupling 
constant. This reduces the sensitivity of observables to the factorization scale fip. Be- 
cause of its importance ZZ production has also been studied to NLO accuracy in the 
SM [i8),j9j]. These results were subsequently updated in ||10llll| . These studies provide 
the precise estimate of higher order effects through K factor as well as the sensitivity 
of the predictions to factorization scale. Importantly, the corrections turned out to be 
larger than the expectations based on soft gluon effects justifying a full-fledged NLO 
computation taking into all the processes. The significance of NLO computations in 
the extra dimension models for the two photon production |12[|13n and Drell-Yan pro- 
duction IfTHITSl has already been demonstrated. Although NLO results are available 
in SM, they do not exist in literature in the context of ADD model for Z boson pair 
production, which is the material of the present paper. 

The results which are presented in this paper are obtained using our NLO Monte 
Carlo code (which is implemented on FORTRAN 77) that can easily accommodate 
any cuts on the final state particles and obtain various kinematical distributions. Our 
code is based on the method of two cutoff phase space slicing (for a review of the 
method see [l6j) to deal with soft and collinear singularities in the real emission con- 
tributions. The idea is to separate the phase space into soft and collinear regions using 
small dimensionless parameters 6s and 5c respectively. These singularities that appear 
as poles in e in dimensional regularization = 4 + e) in the soft and collinear regions 
either cancel with virtual contributions or are mass factorized into parton distribu- 
tions functions. We have evaluated all the soft and collinear sensitive pieces coming 
from real emission and virtual processes analytically. Traces of Dirac gamma matri- 
ces and reduction of one-loop tensor integrals that appear in the intermediate stages 
of the computation are systematically performed using the symbolic manipulation 
program FORM |[T7|. The 75 matrices that appear in the intermediate stages of the 
computation require special care as they are not defined in arbitrary dimensions. We 
have used naive anti-commutation relations between 75 and other gamma matrices 
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in n dimensions and the resulting traces are then computed in n dimensions as they 
are free of 75. Alternatively, one can use other method namely HVBM-scheme which 
was proposed in IITSB and generalized in [19]. In this approach. Gamma matrices and 
momenta in the loop and final state phase space integrals are split into a 4 and an 
n — A dimensional part. The 75 anti-commutes in 4 dimensions and commutes in — 4 
dimensions with rest of the 7 matrices. The results obtained this way coincide with 
those obtained using the naive anti-commutation relations (see |!8l|9l). The finite pieces 
are evaluated using Monte Carlo integration. We will discuss this method briefly in a 
later section which will also serve to introduce some of the required notations. 

The paper is organized as follows. First we very briefly introduce the ADD model 
and set up the notations which we use in this paper. Next we discuss the calculation 
in the framework of slicing method and present the analytical results, and finally we 
present the kinematical distributions. 

2 Next-to-leading order computation 

The fact that electroweak symmetry breaking scale of standard model (SM) cannot be 
made stable against quantum corrections within the framework of SM, leads to excit- 
ing possibility of new physics emerging at TeV scales which stabilize the electroweak 
scale. One of the very fascinating possibilities is the existence of (d) extra spatial di- 
mensions. ADD introduced a model |Illl3 in which the SM fields are localized on a 
3-brane and gravity propagates in all the 4 + dimensions as it is the dynamics of 
spacetime itself. In a simple setting all the d extra dimensions are compactified over 
the same scale R. Although the SM fields are localized on a 3-brane, they can feel the 
extra dimensions through their interaction with the tower of massive Kaluza Klein 
(KK) gravitons. The interaction Lagrangian of SM fields with the KK gravitons is 



given by lEDHZll 
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where k = y/lGw/Mpi and the massive KK gravitons are labeled by a rf-dimensional 
vector of positive integers, n = (ni, n2, . . . , Ud). T^^" denotes the energy momentum 
tensor of the SM. The zero mode corresponds to the usual 4-dimension massless gravi- 
ton. For a given KK level n, there are, one spin-2 state, (n — 1) spin-1 states, and 
n{n — l)/2 spin-0 states, and they are all mass degenerate: 



In this paper we will consider only spin-2 KK states. Let us denote by the fun- 
damental scale in 4 + d dimensions and following [2l| we define the relation among 
the gravitational coupling, the volume of the extra dimensions, and the fundamental 
scale as 



Although the coupling k is Mpi suppressed, the fact that there are a large number of 
Kaluza Klein (KK) modes that couple to the SM fields makes the cumulative effect 
significant. Let us denote the sum of spin-2 KK graviton propagators by X'cfi/ then "Dcff 
times square of the coupling can be written as 



The function /(A/ -y/i) depends on the ultraviolet cutoff A on the KK modes and its 
expression can be found in ||2T|. We will identify A with the fundamental scale Mg. 
The vertex Feynman rules for coupling of SM fields to KK gravitons can be found 
in [[T4ll20ll21| . With all the relevant parameters defined and propagators and vertices 
known we are equipped to carry out the calculation. 

As the gravitons couple to Z bosons, PP — j- ZZ can now also proceed through a 
process where gravitons appear at the propagator level. These new channels makes 
it possible to observe deviations from SM predictions if extra dimensions exist. In the 
following we will consider gravitons only at the propagator level and investigate this 
process at NLO level. 



k'R'^ = 87r(47r)'^/2r(d/2)M7(^+2). 



(3) 
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PP — )• ZZ at NLO has three pieces of computation. A LO piece which is a 2 — )■ 2 
parton level process; second is the 2 — )► 2 order as{= g'^/lQn'^) piece which originates 
from loop corrections; the third and final part originates from real emission process 
where in addition to two Z bosons, a parton is also emitted in the final state. Let us 
take up these three pieces in turn. 

2.1 Leading order 

A leading order parton level process has the generic form 

aipi) + b{p2)-^ Z{ps) + Z{p^). (5) 

In SM this proceeds through quark anti-quark annihilation to two massive Z bosons 
as shown in Fig. 9. The coupling of fermions to Z bosons is 



.eTw 



Y (a - Ca^') , (6) 
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where the coefficients C„ and Ca are defined as 

a = Ti -2Qfsm^ew, Ca = Tl, Tw = . . ^ . ■ (7) 

sm Ow cos 6w 

The Q f and denote the electric charge and the third component of the weak isospin 
of the fermion /, and 9w is the weak mixing angle. We give below the matrix element 
squares summed (averaged) over the final (initial) sate spins, colors and polarizations. 
The SM at LO gives order contribution to the cross sections as given in eq. dH]). 
In addition, two more processes are allowed as the KK gravitons can appear at the 
propagator level, gg — t- G* — t- ZZ and gg ^ G* ^ ZZ, as shown in Fig. [12l These qq 
and gg initiated contributions which are of order are given in eqs. (|9|10|) . 



|M(o)|2 
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qq,sm 
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X (n - 2) - 8mHu{t + u) + {tu - m^) (n - 2) (t^ + m^) 



-2m^(n - 12)tu + 2{n - A)t^u^ 



|M(o)|2 



<l<l,9r 



njsm^ - 16m^(t + -u) + tM(3t^ + 2tu + 
-jsm^ - 24m^(t + u) + tu{7t^ + lOtw + 7u^) 



(9) 



|M(o)|2 



(Ar2 _ 1) 128 



X 



72 



(10) 



-256m^(i + li) + 192m^(t + u)^ - 6Am^(t + u)^ ^s^iAm"^ - s) 

(n- 1) 



n 



188m* - 17r - 226tu - \lu^ + GOm' 



32 f 

+- - 44m^ + 40m^(t + - AOmHu{t + u) + 9tu{t + ti 

(n — 2)^ I 

+m*(-9i^ + 26ixi - du^)} + ^-|692m^ - 13i* - 196^^-^ - 362^^^ 

J n — 2 I 

-imu^ - 13u^ - 5Um^{t + u)- Sm^t"^ + 83tu + u^) 
Next we give the interference of SM qq process with the gravity mediated qq subpro- 
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cess. For convenience we will denote -^g^|*sm^i^,gr + c. by |M(°) p 



qq,int' 



|M(o)|- 
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iiMt QANtu 

Am^{n - l)(i + ii) + Aw?{M - 7)tu{t + u) 



(11) 



X 



+tM|(8 - Zn)t^ - 2{n - 4:)tu + (8 - 3n)u^^ 
-m^l - mu + n{t^ + 22tu + 



here sm and gr repesent contributions from standard model, gravity and interference 
of SM with gravity induced processes respectively, s, t and u are the usual Mandel- 
stam invariants, 'D(s) = Vcs/i and m denotes the mass of Z boson. N denotes the 
number of color and these are exact expressions in n dimensions. A factor of 1/2 has 
been included for identical final state Z bosons. 

The parton level cross sections obtained from the leading order matrix elements 
are independent of factorization scale /Xi?. A LO hadronic cross section obtains its de- 
pendence on fip solely from the parton distribution functions. Due to this sensitivity 
to fip the LO predictions are generally regarded as first approximation. To have a re- 
sult which is less sensitive to hf and which also includes missing higher order pieces 
we need to go beyond the leading order. 

2.2 Next-to-leading order 

In Fig. [lOl the order loop diagrams that appear in SM and in Fig. [13] the diagrams 
that have a graviton propagator are presented. Here we consider only 5 flavors of 
quarks and treat them as massless. These diagrams contribute through their inter- 
ference with the leading order diagrams. In general loop diagrams give ultraviolet 
divergences and infrared divergences when the integration over loop momenta is car- 
ried out. We use dimensional regularization (ra = 4 + e) to regulate these divergences; 
these divergences then appear as poles in e. Note however that owing to the gauge in- 
variance and the fact that the KK gravitons couple to SM energy momentum tensor, a 
conserved quantity, this process is UV finite. Various tensor integrals were reduced to 
scalar integrals following the procedure of Passarino-Veltman The 4-point scalar 
integrals that appear in the gg initiated box diagrams were taken from [23[|24|| . The one 
loop matrix elements are recorded below. The finite pieces of matrix element squares 
denoted by a superscript /in are given in the appendix. 
The SM contribution is found to be 



the interference contributions of SM with the gravity mediated processes are 



qq,int 
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qq,int 
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gg,int 



IM^I 



gg,int 



and the pure gravity contributions are 



T(e) |M(o) 



\,,p. + 4(2C(2)-5)|M(o)| 



(14) 



(15) 



99,9r 



16^ 4 /II 4 



|M(o)|2 



99,9r 



+ -f72C(2) + 70-^- 203 )|M(0)p^^^^^ 



(16) 



where 



^ / N 16 12 

T(6) - -7 + -' 



r 1 + - 



(17) 



r(i + e) y^TTfii 

The theory is renormalized at scale fiR. Cp is the Casimir of the fundamental repre^ 
sentation while Ca is the Casimir of adjoint representation in the color group. 



-1 1 

Cf = , Ca = N, Tf^- 



(18) 



N is the color degree of freedom for quarks and iV^ — 1 for gluons. We can now write 



8 



the order as{jj?j^ contributions coming from virtual diagrams as, 

rfcr"''* = as{fi%)dxidx2f{e,fi%s) 



^ ■ 16 4 /II 



+Cf ciaJ^^^,"'(xi,X2, e) (/g,(xi)4(x2) + xi o X2) 

i 

+Ca da^f'^ixu X2, e)(/,(xi)/,(x2)) 



(19) 



Note the appearance of poles of order 2 in e in the one loop matrix elements. These 
correspond to the configurations which are both soft and coUinear simultaneously. 
These double poles cancel when real emission contributions are included, the remain- 
ing simple poles do not cancel completely and are factorized into the bare parton dis- 
tribution functions at the scale jjp- Several checks ensure the correctness of the matrix 
elements. The Z boson polarization sum —g^^ + k^ky/m? does not give rise to nega- 
tive powers of m. Further, for gluon initiated process the gluon polarization sum is 
^9^lv + {k^rijj + k^n^) / k.n where n is an arbitrary light like vector and the results are in- 
dependent of the vector n. The gauge parameter present in the graviton-gluon-gluon 
vertex |21[ does not appear in the matrix element square; this serves as yet another 
check. Furthermore the SM matrix elements are in agreement with the literature JS]]. 

At NLO we also have to include 2 — )• 3 real emission processes. A generic process 
is of the form 

+ KP2) Z{p,) + Z{p,) + c{p,). (20) 

In Fig. [TT] we show the qq and qg initiated real emission Feynman diagrams which 
appear in SM. In addition, in the ADD model the 2—7-3 diagrams with graviton prop- 
agator are shown in Fig. [I4l Here all the three kinds, qq, qg, gg initiated subprocesses 
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occur. The 2 — )• 3 contributions to cross-section reveal the infrared divergences when 
the integral over the final state particles is carried out. As mentioned above, the sum 
of virtual and real emission cross section is finite after mass factorization is carried out, 
we present very briefly below this in the framework of phase space slicing method. 
For more details we refer to the review ||T6l and our earlier work 1131 . 

Using two small dimensionless slicing parameters 5s and 5c the 2—7-3 phase space 
is divided into soft and collinear regions. The soft is defined as the part of phase space 
where the final state gluon is soft and has an energy less than 5s^/s^2/'^ ^ the center 
of mass frame of incoming partons. In this region the cross section simplifies and we 
have 

da'°f^ ~ asdxidx2f{e, fx% s) + 1^^^ + g in^ 5^^ 



Ca daf^{xi,X2,t)fg{xi)fg{x2) 



(21) 



The region complementary to the soft region is hard region and contains collinear 
singularities. This region is thus further divided into hard collinear region (the region 
of phase space where the final state parton is collinear to one of the initial state parton) 
which contains collinear singularities and hard non-collinear region which is free of 



any singularities. After mass factorization in MS scheme, hard collinear region gives 
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the following contribution to the cross section. 



^^HC+CT ^ as{iiji)dxidx2f{e,iiji,s) 



X 



+2 ( -■ ^ + ^ In ^ ) Ag^g+g llF)fqi{x2, /^f) + ^1 O ^2 
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+d&^^J{Xi,X2,e){ 2 ■ ^fg{Xi,HF)fg{x2,HF) 



+ 2 (^--^ + ^ In Ag^g+g f g {X 1 , ^ f) f g {X 2 , ^ f) ^ 



(22) 



Here terms of order 5s have been dropped and factors of 2 appear because both the 
incoming partons can emit gluons. pi2 — (pi + ^2)^ and the other definitions used in 
the above equation are as follows. 

Ag^g^g^ ^^Pgg{z) = 4C^ 111 5, + ^ , 

dz 



^g^g+g — / ^ Pggi^) 



22 8 

yC^--n/T^ + 8C^ln5,) , 



(23) 



and the function fg^g are defined by 

fqix,HF) = I —fq(-,HF)Pqq{z)+[—fg(-,I^F)Pqgiz), 
fg{x,I^F) = / —fqi-,llF)Pgq{z)+ — fg [- , fiF ) Pggiz) , (24) 



with 



and 



P,j{z) = P,,(z) In ( S^^^ ) + 2i^;.(^). 



P,j{z,e)^Pij{z)+ePl^{z) 



(25) 



(26) 
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Let us now add all the order pieces together; the virtual cross-section rfa'"^* in 
Eq.( [T9l) , the soft piece da^°f^ in Eq. ((2T|) and the mass factorized hard collinear con- 
tribution da^^^^'^ as given in Eq. (|22ll . We see that all the poles in e cancel in the 
sum 

5„ /x^) = da-'-* + da'°f\5s, <5J + 4, /^f)- (27) 

We have made explicit the dependence on the slicing parameters and the factorization 
scale and suppressed other variables. The only order piece, da^~^'"^y{5s, 5c), which 
remains to be included is hard non collinear; it is finite as the integration over 3-body 
phase space here does not include soft and collinear regions. Thus, we need to know 
the phase space only in n = 4 dimensions. It is easy to parameterize the momenta of 
particles in the rest frame of the two final state Z bosons and later lorentz transform 
to the laboratory frame. We can take pi to define the z axis, and using pi and p2 define 
y — z plane. 



Pi = ^i(l,0, 0,l), 

P2 = -£"2^1, 0, sin?/^, cos 



P3 = -^^il,/3xSin92sm9i,/3xCos92sin9i,/3^cos9i), 



P4 = '^^ (1, —Px sin 92 sin 9i, — cos 6*2 sin 9i, —f3x cos 6^1), (28) 



where 



Px = \ll , P34 = bs + Pa)\ (29) 

P34 



and smif) > 0. The momentum p^ is determined by momentum conservation. The 
three body processes are characterized by five independent scalar quantities: 

P12 = (pi + P2)^ = (pi - P5)^ P25 = iP2 - Pbf, 

Pl3 = {Pi - PaY, P24. = {P2 - Pif- (30) 
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Of course any other set of five independent scalars can be chosen. For convenience let 
us introduce the variables x and y, where x = p^a/pii and y is the cosine of the angle 
between pi and P5. We have 

4m^ 



P12 

and 



< X < 1, -l<y<l, (31) 



Pi5 = -^(l-a;)(l-2/), p,, = -P-f{l-x){l + y). (32) 
The phase space is given in 4-dimensions by 

= jt-tt; -^dcos9idx^{l - x)dyd92. (33) 

The integration over the 3-body phase space is carried out using Monte Carlo, and it is 
constrained to avoid collinear and soft regions. The qq and gg initiated processes con- 
tain both kinds of divergences so the integral is constrained using Sg and 4 to avoid 
these regions. The qg initiated process, however, contain only collinear singularities 
(as soft fermions do not give singularities) and the 3-body integration is constrained 
using only 5c. 

We want to express the momenta in the laboratory frame so as to facilitate imple- 
mentation of any cuts, such as rapidity cut, on the final state particles. The transfor- 
mation matrix can be obtained by first boosting to the pi + p2 rest frame and rotating 
to align pi and p2 parallel to the 2;— axis. Finally we will boost to the laboratory frame. 
Successively carrying out these transformations we obtain the final transformation 
matrix Mij whose components are as given below: 

^ X1E2 + X2E1 ^ -2Ei{xiE2 - X2E1) + X2xlS ^^^^ 

XiX2\/S 2EiXiX2\/S 

^ _ Xl^J XiX2S{4:EiE2 - X1X2S) ^ _ y/xiX2S{4:EiE2 - X1X2S) ^^^^ 
2EiXiX2^/S ' X1X2S 

-2Ei(xiE2 + X2E1) + X2xlS X1E2-X2E1 

Mo3 = — 7= , M30 = , (36) 

2EIX1X2V0 XIX2V0 
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M23 = -M20, M32 = Mo2. 



(37) 



The remaining matrix elements are zero. S denotes the center of mass energy of the 
colliding hadrons. The energies Ei and E2 and cos ip can be expressed as 

cos^ = l-f^, (38) 



4V X 



E2 = y^^ 2-{l-x){l + y) . (39) 



Applying the above transformation we can obtain all the momenta in the Laboratory 
frame and can impose any restrictions at the Monte Carlo level. We do not give the 
matrix elements for 2 — 3 processes as the expressions are large and can be obtained 
on request. 

The NLO result is sum da^° +da'^-^°'^y{5s, 5c, fiF)+d(T^-^°'^y{Ss, Sc). The sum da^'^^'^y 
{5s, 5c, lip) +da^~^°'^y{5s, 5c) constitutes QCD correction, but da'^-^°'^y{5s, 5c, Hf) and 
da^~^°'^y{5s, 5c) independently are not physical quantities as these depend on the (arbi- 
trary) slicing parameters. The sum of these two pieces should be independent of the 
slicing parameters as these were introduced at the intermediate stages of calculation. 
A verification of this, in the next section, will serve as a test on the code as to the 
correct implementation of the phase space slicing method. 

3 Results 

In the previous section we have given all the relevant analytical results, now we pro- 
ceed to determine some kinematical distributions. First we demonstrate that the sum 
of 2-body and 3-body contributions is fairly independent of the slicing parameters. 
In Fig. [T] (for SM) and Fig. |2] (for signal) we show the variations of these two pieces 
with the slicing parameters in invariant mass, Q = a/ {ps + PaY, distribution at a value 
of invariant mass equal to SOOGel^. Here both 5s and 5c are varied together with the 
ratio 5s/ 5c fixed at a value of 100 IIT6B . We note that the sum of 2-body and 3-body 
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contributions is fairly stable against variations in these parameters and this gives us 
confidence in our code. In what follows we will use Sg = 10^^ and 5c = 10~^. 

Below we present various distributions for the LHC with a center of mass energy 
of 14 TeV as a default choice. However we will also present some results for a center 
of mass energy of 10 TeV for the LHC. For numerical evaluation, the following SM 
parameters [|25 | 1 are used 

m = 91.1876 GeV, sin^ Ow = 0.231 (40) 

where 9w is the weak mixing angle. For the electromagnetic coupling constant a 
we use a^^ = 128.89. CTEQ6 Il26ll27|| density sets are used for parton distribution 
functions. 2-loop running for the strong coupling constant is used which is given by 
the expression. 



.2 



/3oln 



A 



2 



(/3i//3o) Inln-I^ 

^^QCD 

/3oln- 



QCD I ^'^QCD 



A 



2 



(41) 



with 



M A 

/3i = -Cl--nfTj{?,CF + WA). (42) 

where the symbols are as given in Eq. [18] with number of colors equal to three. The 
number of active light-quark flavors is denoted by n/ (= 5) and the value of Aqcd is 
chosen as prescribed by the CTEQ6 density sets. At leading order, that is at order a°, 
we use CTEQ6L1 density set ( which uses the LO running ) with the corresponding 
^QCD = 165 MeV . At NLO we use CTEQ6M density set ( which uses 2-loop running 
as ) with the Aqcd = 226 MeV; this value of Aqcd enters into the evaluation of the 
2-loop strong coupling. The default choice for the renormalization and factorization 
scale is the identification to the invariant mass of the Z boson pair ie., /i^ = /^i? = <5- 
Furthermore the Z bosons will be constrained to satisfy \yz\ < 2.5, where yz is the 
rapidity of a final state Z boson . 
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In Fig. |3] we have plotted the invariant mass distribution both for the SM and 
the signal, in the range 300 GeV to 1600 GeV. The distribution is presented at the 
higher values of Q as it is in this region the deviations from SM are more pronounced. 
In this plot we display for three extra dimensions ie., d = 3, for fundamental scale 
equal to 2 TeV. We see that for this choice of parameters the signal starts deviating 
significantly from the SM predictions around 400 GeV. To highlight the importance 
of QCD corrections we have also displayed the LO results of SM and the signal, and 
we observe that the K factors (defined as K = da^^'^ /da^^) are large. For the signal 
the K factor is 1.98 at Q = 600 GeV and 1.82 at Q = 1600 GeV. 

To estimate the effect of the number of extra dimension on the invariant mass 
distribution, we plot in Fig. |4] the signal for three different values of d (3,4,5) with Mg 
fixed at 2 TeV. We note that the lower the value of d more is the strength of the signal. 
Next in Fig. |5]we have plotted da/dQ for three different values of Mg (2.0, 2.5, 3.0) at 
a fixed value 3 for the number of extra dimensions. As expected, with increase in the 
fundamental scale the deviations from SM predictions become less, and significant 
deviations from SM are observed at higher energies still. 

If Fig. |6]we have plotted the rapidity distribution da/dY at LO and NLO both for 
SM and the signal for d = 3 and d = A. The rapidity Y is defined as 

Y = -ln^, (43) 
2 P2-g 

where Pi and P2 are incoming proton momenta and g = J53 +J94 ie., sum of the Z boson 
4-momenta. We have plotted this distribution in the interval —2.0 < Y < 2.0 and have 
carried out an integration over the invariant mass interval 900 < Q < 1100 to increase 
the signal over the SM background. As expected the distribution is symmetric about 
Y = 0. 

We have mentioned before that the NLO QCD corrections reduce the sensitivity 
of the cross sections to the factorization scale fip', this we now show in the Fig. [71 We 
have plotted SM and the signal both at LO and NLO, and have varied the factoriza- 
tion scale fip in the range Q/2 < fip < 2Q. The central curve in a given band (shown 
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by the dotted curves) correspond to yu^? = Q. In all these the renormalization scale is 
fixed at fiR = Q. We notice that the factorization scale uncertainties in SM are reduced 
compared to the signal. This is because of the dominant role of the gluon gluon ini- 
tiated process in the signal. Most importantly we have been able to demonstrate that 
a significant reduction in factorization scale uncertainty is achieved by carrying out a 
full NLO computation. For instance at Q = 1400 GeV varying between (5/2 to 2(5 
shows a variation of 20.4% at LO for the signal, however the NLO result at the same 
Q value shows a variation of 6.4%. 

At the end we present in Fig. ^ da / dQ for LHC with a centre of mass energy of 
10 TeV at NLO both for SM and signal. For comparison we have also plotted the 
14 TeV results in the same figure. 

4 Summary and Conclusions 

In this paper we have carried out a full NLO QCD calculation for the production of 
two Z bosons at the LHC at lATeV in the large extra dimension model of ADD. Here 
we take all order contributions, both in the SM and in the gravity mediated pro- 
cesses and their interferences, into account. We have presented all the leading order 
and one loop virtual matrix element squared for the process. The method of two cut- 
off phase space slicing, on which our monte carlo FORTRAN code is based, is very 
briefly discussed. After offering some checks on our monte carlo code we obtained 
invariant mass and rapidity distributions both at LO and NLO. We use CTEQ 6L1 
and CTEQ 6M parton density sets for LO and NLO observables, respectively. Signif- 
icant enhancements over the LO predictions are observed. The K factors are found 
to be large in the invariant mass distribution; for instance the signal has a K factor of 
1.98 at Q = 600 GeV and 1.82 at Q = 1600 GeV. We have also presented the effects 
of variation of number of extra dimensions d and the fundamental scale Mg in the 
Q distribution. We have shown that a significant reduction in LO theoretical uncer- 
tainty, arising from the factorization scale, is achieved by our NLO computation. For 
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instance atQ — 1400 GeV varying fxp between Q/2 to 2Q shows a variation of 20.4% 
at LO for the signal, however the NLO result at the same Q value shows a variation 
of 6.4%. Thus our NLO results are more precise than the LO results and suitable for 
further studies for constraining the parameters of the ADD model. Invariant mass 
distribution is also presented for LHC at a center of mass energy of lOTeV at the NLO 
level. 
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Figure 1: Variation of 2-body and 3-body contributions (of da /dQ at Q = 800 GeV 
SM) and their sum with 6s- Here 6s/ 6c = 100 has been used. 
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Figure 2: Variation of 2-body and 3-body contributions (of da/dQ at Q = SOOGel^ 
signal) and their sum with 5s- Here 5s/ 6^ = 100 has been used. 
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ZZ production at the LHC (V S =14 TeV) 
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Figure 3: Invariant mass distribution at LO and NLO in SM and for the signal at 
Ms = 2TeV and 3 extra dimensions. 
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ZZ production at the LHC (V S =14 TeV) 




Figure 4: Effect of variation of number of extra dimensions in invariant mass distri- 
bution. The fundamental scale Mg has been fixed at 2 TeV. The curves correspond to 
NLO results. 
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ZZ production at the LHC (V S =14 TeV) 




Figure 5: Effect of variation of the fundamental scale in the invariant mass distri- 
bution. The number of extra dimensions has been fixed at 3. The curves correspond 
to NLO results. 
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LHC (V S =14 TeV) 
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d variation, l\/l3=2TeV Y 

Figure 6: Rapidity distribution for Mg = 2TeV for SM and signal for d = 3 and d = A. 
The dotted curves correspond to the LO and solid curves to NLO. We have integrated 
over the invariant mass range 900 < Q < 1100 to enhance the signal. 
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M^P variation, LHC (V S =14 TeV) 




600 700 800 900 1000 1100 1200 1300 1400 
d=3, M3=2TeV Q 

Figure 7: Factorization scale variation in the invariant mass distribution. The number 
of extra dimensions d = 3 and the fundamental scale Mg = 2TeV have been chosen. 
LO curves correspond to CTEQ 6L1 density sets and NLO curves to CTEQ 6M sets. 
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LHC at V S =14 TeV and V S=10 TeV 
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Figure 8: Invariant mass distribution at NLO for SM and the signal. Here the thicker 
curves correspond to = lOTeV and lighter curves to = UTeV at the LHC. 
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5 Appendix 

/^^'^tL f3m^(-4t + u) + 6mH(2t + u) - 2mH'^u(4t + u) 

[t — m^y \ 

+t^u{2t + 3u) + mH{-2t^ + tu- 3w^)) 

+u F2{t) (m^{8t - 2u) - AmH{2t + u)+ t{2t^ + 2tu + u^)^ + t O w| 

+2C(2) - 2Qm^tu{t + u) - Am^if - 8tu + -u^) + ^^(5^^ + Atu + 5u^)^ 

+ , — Sm^fu^it + li) + I2m^it^ - 8tu + u^) 
(4m2 - s)2 V 

+ , ^^f^ ,o f - 32m^tii - 64m^(i + li) + 8mHu(t + uf 
s{Am? - s)2 V 



-(i + uf{?,t^ + 4tu + 3^2) + m^(22t^ + 82t2u + mu^ + 22ii^)) 

(iSm^^it^ - 8tu + M^) 

— s) \ 



{t — m'^){u — m'^){Am^ — s) 
+m^{-9t^ + 131t\ + ISltM^ - 9^3) - 7t\^{t^ + t\ + tu^ + u^) 

-AmHu{At^ + 23t\ + 23tM^ + 4m2) + m^(-9t'^ + 14t\ - mt^u^ 



+lAtu^ - 9u^) + mHu{9t^ + 32t\ + 82t\'^ + 32tu^ + 9m^)) 



(44) 
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g2^2 ^ ^ ^ 



I , ^}' , f - 9m^ - 4m^i2 + 3m^(5i + - m^itifQi + (45) 

+t'^u{2t + 3m)) + (^6m^ - 2m\8t + u) + mH{10t + 7u) 

-t{2t'^ + 2tu + -u^) j + t^uj - - 18m^{t + u) 

-37mHu{t + u) + m^{6f + 80tu + Gu^) + tu{7t^ + Atu + lu^)^ 

G / 
^ 9m^°(i + w) + mHu{t^ + lOiw + w^) 



tu [u — m?) {tv? — t) 
-6m^ {2t^ + btu + 2ii^) + ^^14^(8^2 + Uu + 3ii^) 
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\2 



s {im? — sY 

+ {t + uf{3f + Atu + 3u^) - m^(t + uf{bf + mu + bu^) 

-2m^{lt^ + 31t\ + 37t^/^ + lu^)) + ,,^!^ ( - lUm^it + u) 

^ V itt (4m2 - s) V ^ ^ 

+m^(-19i^ + 334iu - I9v?) - mHu{b3t^ + imu + 53xi^) 

+ ^ 2 ( - 36m^°(t + m) - 8m*^(3t^ - IStw + 3^^) 

tu {i.m? — s) ^ 

+tu{t + uf{3t^ + 4tM + 3m^) + m\3t^ + Glt^w + 61tu^ + 3u^) 
+m^{3t^ - 2fu - mW - 2tv? + 3u^) - 7mHu{t + uf^ 
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where 
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TleVs In ; ^ 



G2 = 47^e^>sL^2 ( At 1 + 27^eX>, In 



In 



{w? - 1) 



+2ImVs7r In 



TZeV., In ( ^ 



2XmV, 



G 



G4 = ZmVsTT, 
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1 62/^2 
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where 
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Figure 9: Leading order diagram in SM. The diagram with the momenta of final state 
Z bosons interchanged (which is not shown here) also contributes. 
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Figure 10: Order virtual diagrams in SM. The diagrams with the momenta of final 
state Z bosons interchanged (which are not shown here) also contribute. 
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Figure 11: Order real emission Feynman diagrams in SM. The diagrams with the 
momenta of final state Z bosons interchanged (which are not shown here) also con- 
tribute. 
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Figure 12: LO gravity mediated diagrams. 




Figure 13: Order gravity mediated virtual correction Feynman diagrams. 
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Figure 14: Gravity mediated real emission diagrams which contribute at NLO. 
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